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Effects of plasma inhomogeneity on self-modulating proton bunches and accelerated electrons 
were studied numerically. The main effect is the change of the wakefield wavelength which results 
in phase shifts and loss of accelerated particles. This effect imposes severe constraints on density 
uniformity in plasma wakefield accelerators driven by long particle bunches. The transverse two 
stream instability that transforms the long bunch into a train of micro-bunches is less sensitive to 
density inhomogeneity than are the accelerated particles. The bunch freely passes through increased 
density regions and interacts with reduced density regions. 
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I. INTRODUCTION 

Studies of plasma wakefield acceleration concentrate 
on various regimes of wave excitation depending on par- 
ticle beam parameters. Early studies of plasma wakefield 
acceleration dealt with linear or weakly nonlinear wakes 
[3, 0] and were supported by proof-of-principle experi- 
ments with low density electron bunches Then 
came the era of the blowout regime stimulated by 
experiments at SLAC [U EH- I n parallel, interest in 
strongly nonlinear positron-driven wakes came up as in- 
tense positron bunches became available [l3l - [l5j |. 

Recently, the concept of proton-driven wakefield accel- 
eration was introduced [3, E3] motivated by the huge 
particle energy and energy content of existing proton 
bunches. It was soon realized [l8j that it is much eas- 
ier to excite the wave by a self-modulating long bunch 
rather than to compress the proton bunch longitudinally 
to sub-millimeter scales. If properly seeded, the trans- 
verse two-stream instability forms the bunch into equidis- 
tant micro-bunches spaced one plasma wavelength apart 
[l9l [2fj| . The micro- bunches then resonantly drive the 
plasma wave to high amplitudes |2l| . 

The process of resonant wave excitation assumes exact 
coincidence of the plasma wavelength and the bunch-to- 
bunch distance. This requirement is automatically met if 
the bunch self-modulates in a perfectly uniform plasma. 
If the plasma density, n, is not constant along the line 
of bunch propagation, then the plasma frequency lo p = 
^jAime 2 /m and the wavelength A p = 2ttc/u! p vary as 
well. Here e is the elementary charge, m is the electron 
mass, and c is the light velocity. It has been shown in 
p3 | that even a moderate density gradient of 0.5% over 
one meter changes the wake phase velocity drastically. 



This is due to the accumulation of the wake phase shifts 
over the long proton bunch. The questions 'how precisely 
must the plasma uniformity be kept?' and 'what are the 
effects of possible non-unformities on the proton bunch 
modulation and electron acceleration?' thus arise. 

We address these questions by two-dimensional simu- 
lations of beam dynamics in plasmas with various density 
perturbations. In Section[II]we define the problem, show 
the overall effect of plasma inhomogeneity on spectra of 
accelerated electrons, and study the effect of non-uniform 
plasmas on the proton bunch. In Section IIIII we look at 
electron dynamics and explain the observed results by 
shifts of the wakefield phase. In Section IIVI we summa- 
rize the main findings. 



II. EFFECT OF DENSITY PERTURBATIONS 
ON THE INSTABILITY 

As a reference point for the study we take the half-cut 
proton bunch (SPS-LHC) discussed in [2l|. This case 
is one of possible scenarios for a proposed experiment 
on proton-driven wakefield acceleration at CERN. The 
initial bunch density is taken in the form 
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with the number of particles N b — 5.75 x 10 10 and sizes 
<j r = 0.2 mm and a z = 12 cm. Other important bunch 
parameters are the proton energy, W b = 450 GeV, and 
the emittance, 0.008 mm mrad. 
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FIG. 1. (Color online) Examples of density perturbations 
tested. The arrows mark the location of the electron injection 
point. 
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FIG. 2. (Color online) Spectra of accelerated electrons for 
density perturbations of various scale (L) and amplitude (Sn). 



The nominal plasma density is no = 7 x 10 14 cm -3 , so 
that c/ujp w ay. The bunch is thus very long if measured 
in plasma wavelengths: oy ps 100 \ v . 

The wakefield is sampled by a side injected low density 
electron bunch. The energy of electrons is 10 MeV, and 
the sizes of the bunch are a e r = 0.2 mm and <r e . 2 = 5 mm. 
The normalized emittance is 2.5mmmrad. The angle 
between the proton and electron bunches is 5 mrad, and 
the beam axes cross at Z = 5.36 m (the proton bunch 
enters the plasma at z = 0). This combination of in- 
jection position and angle was found to provide the best 
quality electron beam at the end of a 10 m long plasma 
column. The electron bunch is timed to get trapped ap- 
proximately 24 cm, or 2ay , behind the leading edge of the 
proton bunch. 

The evolution of the bunches in the plasma is simulated 
by the two-dimensional axisymmetric quasi-static code 
LCODE [23|,[24j with the fluid solver for plasma electrons 
and the kinetic model for beam particles. The applica- 
bility of quasi-static codes to studies of long bunches in 
short scale plasmas is discussed in Ref. po] ]. 

We first perturb the whole plasma column by taking 
its density in the form 



n = 7io[l + 8nsm{2iTz / L)] 



(3) 



(Fig.HJ. The energy spectra of accelerated electrons at 
z = 10 m are shown in Fig. [21 We see that rather small 
perturbations of the plasma density are sufficient to re- 



duce the acceleration rate or even destroy the electron 
bunch. The shorter scale perturbations are more dan- 
gerous than the longer scale perturbations. The picture 
emerging from these spectra is rather complicated and 
results from a mixture of several effects. 

To characterize the effect of density perturbations on 
the proton bunch, we plot in Fig. [3] the maximum wake- 
field amplitude reached at different longitudinal locations 
in the plasma. We see that the instability is rather in- 
sensitive to regular density perturbations. Density varia- 
tions as strong as 2% do not suppress the instability, but 
just reduce the growth rate while keeping the saturated 
wave amplitude at approximately the same level. 

Comparison of the wake amplitude patterns and den- 
sity profiles show that the wake amplitude is largest in 
the regions of depressed plasma density, while in the re- 
gions with increased plasma density, the wake is damped. 
The time dependence of the wakefield amplitude at fixed 
points in the plasma (Fig.@| indicates that the bunch is 
modulated at the plasma frequency of the lower density 
plasma regions. When the modulated bunch enters the 
higher density regions, the modulation period is not res- 
onant and a beating pattern emerges as seen in Fig-Ht 1 - 
This pattern is common to all simulated cases. 

It is instructive to compare the simulation results with 
the theory of instability growth in inhomogeneous plas- 
mas [2g|. Two regimes of the instability are distin- 
guished. For the long scale perturbations with 2ttz/L = 
kjz <C 1, the instability is predicted to be suppressed for 
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where £ = z — ct is the co-moving coordinate measured 
from the bunch head. Our longest perturbations with 
L = 20 marginally fit this regime for z<3m, For this z 
and £ = —24 cm, formula (|4]) gives the cutoff value Sn — 
0.6%. If we take into account that the on-axis bunch 
density is not constant between £ = and £ = —24 cm 
and gradually decreases from n^o to some value, then 
the cutoff perturbation will be even lower. However, in 
simulations we observe almost unimpeded growth of the 
wakefield for Sn — ±2%. In the opposite case of short 
scale inhomogeneity, the cutoff perturbation is 



nbom^z 2 
2n W b e 



1/3 



(5) 



Substituting z = 5m and £ = —24cm into this expression 
yields Sn = 0.8%, which is again well below values (±5%) 
for which the instability is suppressed in simulations. 

The micro-bunches formed in an inhomogeneous 
plasma can be well suited for further wakefield excita- 
tion and electron acceleration. To illustrate this, we take 



3 



+2% 




8n=0 
-0.5% 



> 



L=6.67m -°- 5< 
+0.5% 




+0.5% 




> 



T =0 fim - +20/o 

U - 0m -0.5% flA 
+0.5% . A 




FIG. 3. (Color online) Maximum longitudinal electric field generated by the proton bunch in non-uniform plasmas versus the 
longitudinal coordinate z. Insets under the graphs show corresponding profiles of the plasma density for positive Sn. Crosses 
mark the electron injection point. The circles in the lower right plot indicate the cases shown in Fig. [4] 



the plasma density of the form 
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( 6 ) 

This density perturbation (FigjS^,) contains spacial har- 
monics of different scales in the first half of the plasma 
cell to model the plasma inhomogeneity. In the second 
half, the density is kept uniform to distinguish the effects 
from plasma inhomogeneities on the proton and electron 
bunches. The instability still develops (FigHb), and spec- 
tra of accelerated electrons are almost as good as for the 
uniform case (Fig[5};). This numerical experiment also 
showed the importance of a proper injection point and 
angle. Adjustment of these two parameters was neces- 
sary to provide the best result for the particular instabil- 



ity case. 



III. EFFECT OF DENSITY PERTURBATIONS 
ON ACCELERATED ELECTRONS 

Plasma nonuniformity can affect the quality of the ac- 
celerated electron bunch in three ways: by changing pro- 
ton bunch dynamics and the instability growth, by chang- 
ing trapping conditions, and by changing accelerating 
properties of the wakefield. The instability growth has 
been analyzed in the previous section and was shown to 
be relatively insensitive to density variations. Trapping 
of side injected electrons, even in the uniform plasma, is a 
complicated process deserving a separate study, so we as- 
sume the electron bunch is already trapped and examine 
the acceleration only. We note that the trapping process 
produces a series of microbunches of electrons spaced at 
the plasma wavelength where the trapping has occurred. 
To understand the effect of changing plasma density on 
these microbunches, we perturb the downstream part of 
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FIG. 4. (Color online) Longitudinal electric field in the non- 
uniform plasma with L = 0.6 m and 8n = 2% versus time 
at locations where (a) the wakefield is damped (z = 8.55 m) 
and (b) peaked [z = 9.45 m). Note the difference of vertical 
scales. 



the plasma column: 



no, 
n [l 
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with zq — 6 m. In the first 6 meters of the uniform 
plasma, the self-modulation fully develops and electrons 
are trapped and pre-accelerated to the energy of about 
140 MeV. 

The main effect of density nonuniformity comes from 
dephasing of the wakefield and electron bunches. To il- 
lustrate this, we show in Fig.|5]how the wave phase moves 
with respect to the speed-of-light frame or, equivalently, 
to highly relativistic electrons. Lines in Fig.|6] are coor- 
dinates of the "stable points" where E z = on axis and 
the wakefield is focusing for electrons. Only the wake- 
field fragment near the electron bunches is shown. In the 
first 4-5 meters of the plasma, the phase moves back- 
ward as expected for the developing instability [22|, |27| 
and the wave is unusable for electron acceleration. Af- 
ter the bunch fully self-modulates at z ~ 5 m, the phase 
velocity in the uniform plasma approaches light veloc- 
ity, and the wave can accelerate electrons. In nonuni- 
form plasmas, the phase responds to density variations 
by forward or backward shifts. This can be fatal for 
the accelerated bunches, as illustrated by Fig.[7J In the 
unperturbed plasma, electrons are located somewhere in 
the region of both focusing and acceleration [Fig.[7Ja)]. 
If the plasma density increases, the wavelength shortens, 
the wave near the electron bunch shifts forward, and the 
bunch either enters the stronger accelerating field or falls 
into the defocusing region and gets scattered [Fig.jjfb)]. 
If the plasma density decreases, the wavelength elon- 
gates, and the bunch sees a lower accelerating field or 
even a decelerating field [Fig.[7[c)]. 
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FIG. 5. (Color online) (a) Plasma density perturbation given 
by formula ([6]), (b) longitudinal electric fields generated in the 
uniform and this non-uniform plasmas, and (c) corresponding 
spectra of accelerated electrons at z = 10 m. 



Energy spectra of electrons at z — 10 m are consis- 
tent with this model. For long scale non-uniformities 
[L = 8m, Fig.jHfa)], the density perturbation has a defi- 
nite sign, and we observe either monotonous reduction of 
the beam energy as \6n\ increases (for 8n < 0), or some 
growth of the beam energy followed by complete loss of 
electrons (for 8n > 0). For short scale non- uniformities 
[L = lm, Fig.[5[b)], effects of different signs of Sn are 
averaged, and we observe a simultaneous decrease of the 
electron energy and population. 

To find out the acceptable magnitude of density per- 
turbations <5n max , we refer to the snapshot of the wake- 
field amplitude at the stage of the fully modulated pro- 
ton bunch (Fig. [5]). There are two important points on 
this graph. One is the location of the electron bunch 
(£ ~ —24 cm). Another is the place in which the wake- 
field is mainly formed, or more strictly the point where 
the field has a strength of half of the field that accelerates 
electrons (£ ~ — a z = —12 cm). The distance between 
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FIG. 6. (Color online) Locations of constant phase points (where E z — on axis and dE z /d£, > 0) versus propagation distance 
for the uniform (thick line) and perturbed (thin lines) plasmas. Grey bars show locations of trapped electrons. Insets under 
the graphs show corresponding profiles of the plasma density for positive Sn. 
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FIG. 7. (Color online) Phasing of the accelerated electron bunch in plasmas of the proper density (a), increased density (b), 
and reduced density (c). 
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FIG. 8. (Color online) Spectra of accelerated electrons for the 
perturbed downstream part of the plasma column. 
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FIG. 9. (Color online) Instant snapshot of the wakefield am- 
plitude for ct = 5 m. 



the two is the "length" of the useful wakefield. Let N pel 
be this distance measured in plasma wavelengths. In our 
case iVper ~ 100. Since \ p oc n -1 / 2 , then the relative 



change 5n of the plasma density results in the relative 
change SX p — —Sn/2 of the plasma wavelength and shift 
of the wakefield phase by the fraction N pcl Sn/2 of the 
period. If the electron bunch is located in the middle 
of the focusing and acceleration interval, as is shown in 
Fig.JTJa), then the shift of the wave by 1/8 of the pe- 
riod will bring the bunch to a defocusing or decelerating 
phase. We consider this shift as the maximum acceptable 
one. Thus we come to the engineering formula 

5n max ~ 0.25/iV pcr . (8) 

Here <5n max ~ 0.25% in agreement with Fig.[5J 

If the density perturbation is much stronger than the 
acceptable value (J5J), then the electron bunch is scattered 
away at the defocusing distance 



(9) 



where j e is the relativistic factor of the electrons, and E 
is the wave amplitude in units of the wavebreaking field 
mcujp/e. The distance © is also the maximum scale of 
the increased plasma density region which the electron 
bunch can pass through. In the case considered, we have 
E « 0.1 at z = 6 m. Thus lOMeV electrons can sur- 
vive in plasma nonuniformities up to 0.25 cm long. For 
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FIG. 10. (Color online) Fraction of the electron bunch remain- 
ing in the circle of radius ay for strong plasma nonuniformities 
of different length scales L. 



140 MeV electrons this length is 1 cm, and for 1 GeV elec- 
trons 2.5 cm. The distances are quite short, so keeping 
the plasma uniformity better than is a must once 
electrons are injected into the wakefield. Figure [TU] il- 
lustrates how quickly electrons are lost for plasma per- 
turbations (J7J with 5n = 0.01 > <5n max and various L. 
For L > Ld, most of the bunch is scattered by the first 
density maxima. For L < Lj, the bunch survives over 
much longer distance under the action of some averaged 
focusing force. 



IV. SUMMARY 

Plasma inhomogeneities result in varying frequencies 
of the plasma oscillations, causing phase differences be- 
tween the driving force (the modulated proton bunch) 
and the plasma wave, as well as phase differences be- 
tween the accelerated particles and the wakefields. If 
accelerated particles are located N pcl - wakefield periods 
behind the location of wakefield formation, then relative 
density perturbations stronger than ([8]) are sufficient to 
destroy the accelerated bunch at the defocusing distance 
(|9]). This effect imposes severe constraints on density uni- 
formity in plasma wakefield accelerators driven by long 
particle bunches. 

The transverse two stream instability that transforms 
the long drive bunch into a train of micro-bunches is 
less sensitive to density inhomogeneity than are accel- 
erated particles. The drive bunch freely passes through 
increased density regions and interacts with reduced den- 
sity regions. However, the spacing of the micro-bunches 
produced in this case corresponds to a plasma density 
which is lower than the path average. 
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